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We study a mixture of Ni, bosons with point hard-core boson-boson interactions and Nf non- 
interacting spinless fermions with point hard-core boson-fermion interactions in ID optical lattice 
with external harmonic confine potential. Using an extended Jordan- Winger transformation (JWT) 
which maps the hard-core Bose-Fermi mixture into two component noninteracting spinless fermions 
with hard-core interactions between them, we get the ground states of the system. Then we deter- 
mine in details the one particle density matrix, density profile, momentum distribution, the natural 
orbitals and their occupations based on the constructed ground state wavefunctions. We also discuss 
the ground state properties of the system with large but finite interactions which lead to the lift of 
ground degeneracy. Our results show that, although the total density profile is almost not affected, 
the distributions for bosons and fermions strongly depend on the relative strengthes of boson-boson 
interactions and boson-fermion interactions. 

PACS numbers: 03.75.Hh, 67.85.Pq 



I. INTRODUCTION 

Recently, strongly interacting gases of bosons in one 
dimension have been experimentally realized P, [2, Q 
by loading a Bose-Einstein condensate into a deep 
two dimensional optical lattice to create arrays of one- 
dimensional (ID) atomic systems. The achievement of 
the Tonks-Girardeau (TG) regime [3] in an optical lattice 
has stimulated extensive theoretical interest in the study 
of many-body physics of ID quantum gases [1, 0, 0j 9- 
With the increase in the interaction strength, the ID 
Bose gas evolves from a Bose-Einstein condensate to 
"fermionized" TG [1, HH, E El, Q El • The mi- 
croscopic mechanisms of the evolution has been recently 
studied by both analytical methods [H, EH and various 
numerical methods [l^, El, E3, EB| ■ In the limit of the 
infinitely repulsive interaction, the many-body state of 
a TG gas has been shown to correspond to the states 
of a noninteracting Fermi gas via a Bose-Fermi mapping 
For the lattice model of TG gas, the Hamiltonian 
in a periodic lattice can be mapped onto the ID XY 
model of Lieb, Schulz, and Mattis, which has been exten- 
sively studied in the literature [16|. With an additional 
confining potential, the lattice TG gas has been stud- 
ied by means of an exact numerical approach by Rigol 
and Muramatsu [11 El, HI- Very recently, Girardeau's 
Bose-Fermi mapping method has also been generalized 
to deal with mixtures of multi-component quantum gases 
[mils, 23, 24]. 

On the other hand, mixtures of bosonic and fermionic 
atoms have been studied extensively as they initially pro- 
vided a convenient way to achieve degenerate fermionic 
gas by means of sympathetic cooling [H, Due 
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to their rich phase diagram, the Bose-Fermi mixtures 
have attracted many theoretical studi es [27l [281 [29l . [30l . 

im, m [n H m m, 113, m, m, [loirwiif- 

Among those works, particular attention has been paid 
to the ID model of mixed bosons a nd p olarized fermions 
[23, m, m, [13, [m m m, S, [ai, in which there 
are only s-wave scattering for boson-boson and boson- 
fermion interactions present. While most of these inves- 
tigations relied on the mean-field approximations or the 
Luttinger liquid theory, there is rarely exact result ex- 
cept for the homogenous model with equal boson-boson 
and boson-fermion interactions, which is exa ctly solved 
by the Bethe-ansatz method [11, [H, [H, [H, H^. How- 
ever, for a lattice system, there is no exact result even 
for the lattice correspondence of the integrable contin- 
uum Bose-Fermi system (4^ . In this work we mainly 
study the lattice Bose-Fermi gas in the hard-core limit 
where both the boson-boson and boson-fermion interac- 
tions are infinitely strong. In this limit, wc can apply 
an extended Jordan- Wigner transformation and an exact 
numerical approach, which can be viewed as an extension 
of the method by Rigol and Muramatsu [ll|, to treat 
the hard-core Bose-Fermi mixture exactly. We focus on 
the ground state properties and analyze the behaviors of 
the one-particle correlations, the momentum distribution 
function, the natural orbitals and their occupations. The 
properties of the system with large but finite interactions 
are also discussed. Since the interaction can be tuned in 
principle within a large range of regime by exploiting 
Fesbach resonance, our results might be experimentally 
relevant. 

The content of the paper is as follows. In the next 
section we describe the model system first and then the 
approach used to calculate the ground state one parti- 
cle density matrix exactly. In Sec. Ill, we discuss the 
properties of hard-core Bose-Fermi mixture confined in 
harmonic traps. In the Sec. IV, we discuss the proper- 
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ties of the system with large but finite interactions. The 
paper is concluded in Sec. V. 



II. SYSTEM AND METHOD 

We consider a mixture of Ni, bosons with point hard- 
core boson-boson interactions and Nf noninteracting 
spinless fermions with point hard-core boson-fermion in- 
teractions and assume that the boson and fermion par- 
ticles have the same masses mb — mf, which could be 
realized by choosing an isotope of a given alkali element. 
Let Xb = {xib, ■■■,XNi,b) and Xf = {xif,...,XNjf) indi- 
cate the boson and fermion coordinates respectively. For 
the system trapped in the potentials including the optical 
lattice and an additional harmonic trap, the Hamiltonian 
of the system is 

H = Hb + Hp + Hbb + Hbf (1) 



with the from of 

L-l 

HHub = -t^(6l6,+i+///,+i+i/.c.) 
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+^Y.^M-i) + Ubij2<4 (3) 
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where b\ (fj) and bi (fi) denote the bosonic (fermionic) 
creation and annihilation operators at site i, respectively, 
and they satisfy the standard (anti-) commutation rela- 
tions, i.e., 6]] = S^j, {fi, /]} = Sij, and [b.„ fj] 0. In 
the hard-core limit, the Hamiltonian is simplified to 

H = Hb + Hf (4) 
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where Qbet 9bf ~^ °° under the hard-core condi- 
tion which means that the many-body wave function 
"ifiXbjXf) of the system vanishes at all boson-boson 
(BB) and boson-fermion (BF) collision points. Explic- 
itly, the potentials take the form 
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-t ^{blh+i + H.c.) + Va^ J2 
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where additional on-site constraint 

b\h + ///. < 1 

is assigned to avoid double or higher occupancy 'iB^I . Here 
t is the hopping parameter to be decided by the optical 
lattice; L is the number of the sites; V{i) = Va?i^ is the 
harmonic confined potential with a is the lattice space 
and V is the strength; n\ = b\bi and n{ = f} fi are 
the boson and fermion particle number operators respec- 
tively. Given the local Hilbert space at i-th site composed 
of a set of {|0), 6j|0), //|0)} under the single occupied on- 
site constraint, the on-site constraint can be written in 
the follow forms 
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""bjix) = Vbj SIT? [ttx / a) + -mbjuoljX 



where a is the lattice spacing associated with wave vector 
fci = 7r/a of the standing laser light. In this work, we 
consider only the case with the trap acting on bosons and 
fermions being the same, ie., Vb{x) — Vf{x) — v{x) with 



and LUb — ujf. 



We use the Wannier function (only consider the lowest 
Bloch band) of the optical lattice to expand \E'(X{,,Xj) 
and get the second quantized Hamiltonian of H, which is 
the standard Hubbard model of Bose- Fermi mixture 13111 



{b,M} = i-fjf 



1 - bjh, 



{/.,//} 

and the following equations are also valid 



6f 



blfl = fA = hbl = hf} = 0. 



(6) 



(7) 



In order to get the ground state properties of the sys- 
tem, we extend the general Jordan- Wigner transforma- 
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tion (JWT)[l3 and get the following transformations: 
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n 
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which map the Hamiltonian of the bosons into noninter- 
acting spinless fermions Hamiltonian. Using the JWT we 
can change the Hamiltonian of the system into 



with 
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where a =], j, and 
for i ^ j, else 

t r t 1 , , 



(10) 



for the on-site constraints. Here = c\^Cia is the a- 
kind (we note the spinless fermions from the boson by J,- 
kind and the original fermions by |-kind) fermion number 
operator, and iV|(|) = ^fc(/)- The Hamiltonian Hi de- 
scribes a mixture of two component fermions with point 
hard-core interactions between two kinds. Notice that 
the operators anticommute between two kinds. 

Next we construct the ground state of the Hamilto- 
nian Hi under the constraints Eq. fTU)) with the method 
proposed by Batista et al ^7[. We consider a set of par- 
ent states, labeled by the string configuration cr, with 
iV = iV| -t- A^l particles and L ~ N holes, \<^o{d)), and 
the form is: 



l*o('?)) 



cricr2f3 
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where Ui indicates the kind (| or [) of the fermion par- 
ticle at site i, L is the number of sites. Notice that the 
number of the configuration a is C;;\ Then we rewrite 
the Hamiltonian Hi with Hi = T + Hy , and 



T = 



^ZfT ; — C^^Cj-|_l^cr H~ H .C. 



Hv = Va^^fnl 
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The states |$o((t)) are eigenstates of Hy and they are 
degenerate with different a. 

By applying the hopping operator Tia we can generate 
a subspace M{a) from the parent state |$o(f?))) and we 
denote 



or, in general 



|$i('?)) = rjv,,|$o('?)) 



|$,(<7))=T,.|$,(a)). 



Obviously the dimension of the subspace M{(t) is Cj; , 
and there are C^*" subspaces. Moreover these different 
subspaces for different a are orthogonal. 

Next we construct the ground state in the subspace 
A/((t). For a specific a, we can make the following map- 
ping: 
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which maps the two component fermions (cio-) into a sin- 
gle spinless fermion (ci). It is straightforward to show 
that in the corresponding new basis the system Hamilto- 
nian can be written as 
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The ground state properties of the fermionic system 
H spinless with TV particlcs have been analyzed in Ref . [Tsl . 
Following the approach therein, we let P denote the low- 
est N eigenfunctions of the Hamiltonian Hspiniess which 
can be obtained by diagonalizing Hspmiess ■ 



P = 



( Pii Pi2 ■■■ Pin \ 
P21 P22 ■ ■ ■ P2N 
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where Pin are the coefficients of n-th single particle state 
IV'n) = X^i^i -Pmc||0). Then the ground state of the spin- 
less fermion gas is the state with the lowest N eigenstates 
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fully filled, and the form is: 
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= X]det(P«)ct|0) (16) 



where s index the combination formed by taking A'' num- 
bers from the set A = {1,---,L}, P^isa square matrix 
with N ranks that the N rows are taken from P according 
to the combination s. represents cl_^cl^ ■•■cl^, and Si is 
the i-th number in the combination s. We had assumed 



that the numbers in combination are all sorted ascend- 
ing. Now, we can use the reverse mapping of Eq. (jl3p to 
the get the ground state of the Hamiltonian Hi in the 
subspace M{(j) 

\'^hM)) = Edet(A)(-ir'4,^cj,^|0>. (17) 

s=l 

where q index the combination formed by taking Nb num- 
bers from the set T = {1, • ■ -^N} which means that 
the qi{i = 1, • • •,iVb)-th site in a is occupied by j-kind 
fermion and qi is the i-th number in the combination 
q, q is just another way to index a. cj^^ represents 

ct ■ ■ ■ c\ ,q represents the combination T — q, and 

Tq notes the times of the permutation to put the set 
{sqi , • ■ ■, s,„^ , , • • ■, Sq^^ } iuto s. Reminding that the 
forms of Hamiltonian for f and |-kind fermion are the 
same, the ground states I'^H^^i'^)) are C^*" degree degen- 
erate because of C^'' different a. 

Supposing the ground state of Hi given by 1$^''^), the 
one-particle density matrix function of boson of the sys- 
tem can be written in the form: 



pf, = (<fg"^|6l6,|<i>g") 



,BF\ 
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where \^q^) is the ground state wave function of Bose- 
Fermi mixture and 
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In order to calculate ($^^1 (and |<i>^)), it is convenient to 
use the following identities: 



i-l 



and 



/3=1 /3=1 



(3=1 13=1 



(20) 



(21) 



where 2; = 1 if j < i, otherwise z = 0. Following the same 
way shown above, we can get the one particle density 
matrix function of fermion [pf^ = {^Q^\fjfj\^Q^)) and 
other quantities such as correlation functions. 
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FIG. 1: (Color online) Top panels: contour plots of the 
one particle density matrices. p^i^{Left panel), p^^^: (middle 
panel) are the one particle density matrices for the Bose Fermi 
mixture with A^i, = 3 and Nf = 2. Right panel, the one 
particle density matrix(p^„re) for the pure Fermi gas with 
Nf — 5. Bottom panels: corresponding off-diagonal sections 
along the anti-diagonal. All the systems are with 13 sites and 
Va^ = 0.02t, where V is the strength of the harmonic trap 
and a is the lattice spacing. 



III. HARD-CORE BOSON-FERMION MIXTURE 
IN THE HARMONIC CONFINE POTENTIAL 

Since the ground state in the hard-core limit has a 
degeneracy of C^'' , for convenience we first consider the 
case that \^q^) is formed by the summation of l^'^j (q)) 
with all the degenerate states having the same weight, 
i.e.. 



i==EEdct(^^)(-i)'''i/Uioi22) 



The ground state \^q^) is related to 1$^''^) by the gen- 
eralized JWT. We note that the above construction is 
essentially equivalent to the construction of generalized 
Bose-Fermi mapping by Girardeau et. al 21]. Then fol- 
lowing the method proposed in the above section we can 
work out the density matrix function and show them in 
FiglU which arc found to fulfil the relation 



B 
rmix 



Nb 
N 



P™(*,J) 



(23) 



according to the data, where p^^^ is the one particle 
density matrix of boson for the Bose Fermi mixture and 
P^'^ihj) is the one particle density matrix of pure TG 
gas of N bosons obtained by the method proposed by 
Rigol and Muramatsu ,1^,]. It is easy to see that the one 
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FIG. 2; (Color online) The density profiles for several systems, 
n^'^' {{) is bosonic(fermionic) density profile for a Bose Fermi 
mixture with A'^;, = 3 and Nf — 2. Pa;B stand for the pure 
TG Bose gas with Nt — x, and Pa;F stand for the pure free 
Fermi gas with Nf = x. Again all the systems are with 13 
sites and Va'^ = 0.02t. 



FIG. 3: (Color online) The momentum distributions for sev- 
eral systems. The systems are defined by the same way de- 
scribed in Fig[2] Notice that the momentum k is in units of 
kL which is the wave vector of the optical lattice. 



particle density matrix of boson decays as the distance 
grows, while the fermionic one decays as the distance 
grows but exhibits typical sign changes due to the Fermi- 
Dirac statistics. Same result of Eq. (p3)l for the continuum 
systems have been found by Girardeau et. al 21]. 

The bosonic and fermionic one particle density distri- 
butions n^„(i) = and n^„(i) = p^i^{i,i) are 
both proportional to the density n^'^{i) = p^'-^{i,i) of a 
TG gas of N bosons [1], i.e., 

Nf Nb N N ' ^ ' 

where n^{i) is the density of the noninteracting gas of N 
fermions in trap which is same to 71-^*^(1) [l^. Actually 
under the state the rate of the probabilities that 

the site i occupied by the "f-kind and j-kind fermions is 
Nh/Nf. This result means that there is no phase separa- 
tion between bosons and fermions. We show the numer- 
ical results of the density profiles in Fig[2] which a gree 
with theoretic results obtained by Girardeau et. al [2l[ 
and Fang et. al [13] ■ Comparing distributions of the TG 
gas of 3 bosons and the free fermion gas of 2 fermions 
with the ones for the mixture(n^^^^), we can see that 
as the other kind particles adding in, the origin particles 
have to hold the higher energy states and the density 
distributions become boarder with lower weight. 

The momentum distributions are defined by the 
Fourier transforms with respect to i^j of the one particle 




FIG. 4: (Color online)The occupation of the natural orbitals 
for several systems defined by the same way described in FigH 

density matrices with the form of 

+00 

n^(^)(fc) = |^'(fc)|2 e-*'=("-"V-^^^'(",™), 

71.171 — — 00 

(25) 

where ^'(fc) is the Fourier transform of the Wannier func- 
tion, and k denotes momentum. Since the bosonic one 
particle density matrix for the mixture is proportional to 
the TG one, the bosonic momentum distribution for the 
mixture is also proportional to the TG one. The numer- 
ical results of the momentum distributions are shown in 
FigEl The peak structure in the momentum distribu- 
tion of boson reflects the bosonic nature of the particle, 
and is in contrast with the structure of the momentum 
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distribution for the equivalent noninteracting fermions. 
Again the distributions for the mixture are boarder and 
fower than the pure ones because of holding higher energy 
states. 

The natural orbitals {<l)^'^^\i)) are defined as the 
eigenfunctions of the one particle density matrix [46l|: 

M 

5:P^(^)(*,,)<(^)0-)=A^(^)<(^)(z) (26) 

and it can be understood as being effective one particle 
states with occupation A,,. In FigU we show the oc- 
cupations for boson and fermion. The occupations are 
plotted as a function of the orbital numbers 77, and they 
are ordered by starting from the highest one. As the one 
particle density matrix of boson of the mixture is propor- 
tional to the TG one, the natural orbitals of boson of the 
mixture arc the same as the TG one with the occupations 
fulfilling the following relations 



V 



(27) 



where A™ denotes the occupation of 77-th natural orbital 
for a pure TG gas composed of N hard-core bosons. The 
peak on the lowest orbital of boson is the feature of the 
boson. As for the fermion, the occupation is no longer 
the step function with lowest N orbitals fully filled as the 
distribution of the pure noninteracting A'^ fermions. And 
there is no peak at the lowest orbital. 



IV. THE PROPERTY OF THE SYSTEM WITH 
LARGE BUT FINITE BB AND BF REPULSION 

As we have discussed in section II, the ground state in 
the hard-core limit has a huge degeneracy. However, we 
expect that the degenerate ground state would be lifted 
when the on-site interactions deviate the infinite limit. 
Next we consider the case with Ubb and Ubf being large 
but finite, for which the second quantized Hamiltonian of 
H is the standard Hubbard model of Bose- Fermi mixture 
with the form of Eq.(I3]). In the situation that Ubb and 
Ubf are still large, the state with double occupancy on 
the same site is a high-energy state and we can use the 
standard projection method to derive the low-energy ef- 
fective Hamiltonian of the system which is given by 



L-1 



i=l 

L L 
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1=1 

2 L-l 
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Ubb^ 



(28) 



with Sj = blf,, S, = {n\-ni)/2. In the limit Ubb-, Ubf ^ 
cxD, the last two summation terms in the Hamiltonian 
vanish and the system reduces back to Eq. ([24|) which we 
studied in the last section. As Ubb and Ubf become finite 
but large, the last two summation terms can be viewed 
as perturbations to the system, so one can expect that 
they wouldn't cause significant changes to some prop- 
erties of the total system such as ground state energy 
and total density profile (see Fig[6|3) because of terms 
of /Ubf and /Ubb being very small for large Ubb and 
Ubf- Nevertheless, a significant effect induced by these 
small terms is the lift of the degeneracy of the ground 
states and the true ground state would be a recombina- 
tion of C^' degenerate states with the weights of states 
to be determined by minimization of the energy due to 
perturbation terms. 

One observes that terms of /Ubf lead to an effec- 
tive isotropic antiferromagnetic exchange interactions be- 
tween "spins" (bosons or fermions) on neighboring sites 
to lower the ground state energy. On the other hand, 
terms of /Ubb produce an effective attractive interac- 
tions between neighboring bosons and thus states with 
all the bosons concentrated together have lower energy. 
Consequently, the relative distribution for the bosons or 
fermions will be changed and determined by terms of 

/Ubb and /Ubf to further lower the ground state en- 
ergy. However, one can expect that the total distribu- 
tion shall not be changed too much because the terms 
of /Ubb and /Ubf are very small in comparison with 
hopping terms. As the system Ea. ((28|) has no analyt- 
ical results any more when the interaction parameters 
Ubb,Ubf being finite, as an approximation, we can treat 
the charge part and spin part separately. Then the charge 
part, which does not distinguish bosons or fermions, is 
determined by the matrix P according to Eg. psp with 
the state given by Eg. pop . whereas the spin part, which 
decides the weight of states with different spin configu- 
rations, is determined by Hamiltonian with only Ubf and 
Ubb terms in Ea. l|^5|) . Considering that, with the presence 
of the harmonic trap, the particles mainly concentrate at 
the center of the trap with holes around, the spin part is 
defined on the system with Nb bosons and Nf fermions 
on Nb + Nf sites. This approximation is similar to the 
spin charge separation approximation. 

First, we consider the situation Ubb ^ Ubf ^ 1- In this 
limit, the terms of Ubf dominate and we can set Ubb — 00. 
Because of the spin fluctuation terms of SfSj + SfSj, 
the off-diagonal terms appear between different subspace 
M{q), and the ground state of the system becomes com- 
plicated. Actually in the limit Ubb = 00, the Hamiltonian 
Huub can be mapped to the Fermi Hubbard model by a 
JWT similar to Eq.® and they have the same thermo- 
dynamic properties _48]. As for the ground state prop- 
erties, it can be worked out by the mapping from the 
fermi Hubbard model. When Ubb is away from the in- 
finite limit, there are no analytical results. In Fig|5]we 
show the density profiles in the situation Ubb ^ Ubf ^ 1 
with spin charge separation approximation. From the 
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FIG. 5: (Color online) The density profiles of the Bose Fermi 
mixture with 3 boson and 2 fermion on 11 sites, the trap 
strength Va^ = 0.2t, U^b = 2000t and U^s = 200f. 'ED' 
stands for exact diagonalization. 'EDS' stands of the results 
gained by spin charge separation approximation with the spin 
part (weights of subspace ground states) is determined by 
exact diagonalization. 
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data, the total distribution is the same with that of the 
pure TG boson gas. Because of the effective antiferro- 
magnetic exchanges between bosons and fermions, the 
specie-dependent distributions for bosons and fermions 
exhibit quite different behavior with alternating peaks. 
To check the validity of this approximation method, in 
FiglS] we also show some results of small system of Bose- 
Fermi Hubbard model gained by exact diagonalization 
method [38!|. We can see that the results gained by the 
approximation agree well with the results obtained by ex- 
act diagonalization and the density distributions of the 
total particle are almost the same for the data obtained 
from both methods (the difference < 10~^). 



Next we consider the situation Ubf ^ Ubb ^ 1- Then 
the Ubb terms dominate and for simplicity we can first 
set Ubf — oo. Although the Ubb term is diagonal in 
the Hilbert space UqM{q), the states |*g^(g)) (Eq. ^TT]) ) 
with different spin configurations are not degenerate any 
more for different q because of the Ubb term. Since the 
Ubb terms tend to make the boson concentrated together, 
the ground states of the configuration q with all the bo- 
son staying together have the lowest energy. Then the 



ground states have the degeneracy of 



Nf + l 



^ Nf + 1. 



We suppose that the ground state \^q^) is formed by 
these degenerate states in subspace M{q) with all the 



FIG. 6: (Color online) (a) :The density profiles of the Bose 
Fermi mixture with 3 boson and 2 fermion on 11 sites, the 
trap strength Va'^ = 0.2t, Utb = 200t and Ubf = 2000t. The 
curve with symbol square (circle) is the density profile of the 
boson (fermion) for the Bose Fermi mixture calculated from 
the ground state \^a^) in Eq[29l (b): The total density pro- 
files of difi'erent Ubb, Ubf for the Bose Fermi mixture with 3 
boson and 2 fermion on 11 sites, the trap strength Va^ = 0.2t. 
The curve with symbol square is the total density profile cal- 
culated from the ground state in Ea l22l 



degenerate states having the same weight, say: 

Nf + l 

1*^0 - ^7=^ E (29) 

Nf + l 

= -mm ^ Edet(p.)(-ir"cL,,ct.^io). 

V / 9' = 1 ^=1 

Then we can get the density matrices and other quanti- 
ties. In FigH^, we show the density of the system form 
the state 1$^''^). We can see that there is a phase sepa- 
ration in the system with bosons are in the middle of the 
trap and fermions surround them. As Ubf is away from 
the limit cx), the Nf + 1 fold degeneracy of the ground 
state is split. For comparison, t he density profiles in the 
limit Ubf ^ Ubb S> 1 with spin charge separation ap- 
proximation and exact diagonalization are also shown in 
FiglS^- It is clear that the results obtained by spin charge 
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separation approximation agree well with the ED ones, 
however the results gained by the state \^Q^)(Eq.^^) 
do not agree very well with the others because the Ubf 
terms are neglected. But the results gained by the state 
$^''^UEq. ([29|) ) and spin charge separation approxima- 
tion both indicate that there is a phase separation in the 
system with bosons located in the middle of the trap and 
fermions surrounded. Again the density distributions of 
the total particle are almost the same from the data for 
the three methods (the difference < 10~^). 



matrix, then we yield the density profiles, momentum 
distribution, the natural orbitals and its occupations. We 
also discuss the property of the system with large but 
finite interactions. We find that, despite the total density 
distribution not sensitive to relative strengths of Ubf and 
Ubb, the boson and fermion distributions rely on Ubb 3> 
Ubf » 1 or Ubf » Ubb ^ 1- We hope that our study 
could be helpful for the experimental achievement of the 
ultracold boson- fermion mixtures with hard-core BB, BF 
interactions in optical lattices. 



V. SUMMARY 

In conclusions, we have studied in detail the ground 
state properties of the mixture of the hard-core bosons 
and noninteracting fermions with point hard-core boson- 
boson and boson-fermion interactions in the ID optical 
lattice with harmonic confine potential. Using extended 
Jordan- Wigner transformations, we calculate the density 
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